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A parametric design analysis was performed using a heat pipe in an 
existing deep sea Radioisotopic Thermoelectric Generator (SNAP-21). Heat 
is transferred from an annular fuel pellet to an annular thermoelectric 
generator through a connecting heat pipe. The fuel pellet is fully 
shielded so that the thermoelectric generator is easily removable. 

Overall efficiency and the weight of major components were determined 
for varying fuel radii of from 1.3 inches to 1.7 inches and for varying 
insulation thicknesses of from 1.0 inch to 2.0 inch. 

The analysis indicates that there is a particular fuel radius (at 
constant insulation thickness) at which minimum weight is reached, 
while the maximum overall efficiency is obtained at a larger fuel radius. 
The median design has an overall efficiency (at the beginning of life) 
of 5.4% and a total weight of 570 lbs. These design results, when com- 
pared to the existing SNAP-21 design gives an increase in overall ef- 
ficiency of at least 7%, and a reduction in total weight of 127o. 
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NOMENCLATURE 



Symbol 




B 


nuclear radiation buildup factor 


D 


dose rate due to nuclear radiation (mr/hr) 


E 


energy of nuclear radiation (ev) 


k 

c 


thermal conductivity of the fuel capsule (BTU/hr-f t-°F) 


k f 


thermal conductivity of the fuel (BTU/hr-f t -°F) 


k I 


thermal conductivity of the insulation (BTU/hr-ft-°F) 


k 

s 


thermal conductivity of the biological shield (BTU/hr-f t-°F) 


LF 


length of fuel (in) 


p 


source strength for volume source (Curies) 




heat transfer rate to heat pipe from fuel/unit 
length (BTU/hr-f t) 


( a } out 


heat transfer rate through insulation from fuel/unit 
length (BTU/hr-f t) 


RF 


outer radius of the fuel (in) 


RI 


outer radius of the insulation (in) 


RS 


outer radius of the shield (in) 


t 

c 


self absorption thickness (in) 


t. 

1 


shield thickness (in) 


T 

H 


hot end temperature of thermoelements (°F) 


Z 


distance from source to where dose due to nuclear 
radiation is to be measured (in) 


V 


Carnot efficiency (%) 


^Om 


material efficiency (7 0 ) 


V 


overall efficiency (%) 


yUi 


attenuation coefficient of shields (1/in) 


JJq 


attenuation coefficient of source (1/in) 
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I. INTRODUCTION 



Future demands for measurements of oceanographic data needed to 
support scientific, commercial, and military operations necessitates 
the development of reliable power sources for deep sea use. For these 
undersea missions that require continuous, relatively low levels of 
electrical power for a few years, radioisotope powered systems are 
needed. Because isotope fueled themoe lec trie generators do not require 
oxygen and have no continuously discardable waste (i.e., no exhaust), 
the performance of these systems are independent of the environment in 
which they operate [14]^* These advantages make them the ideal source 
of power for this application. 

The 3-M Company has designed a radioisotopic thermoelectric genera- 
tor (RTG) for just such a purpose (SNAP-21). The basic features of the 
SNAP-21 design are shown in figure 1, with the SNAP-21 system capabili- 
ties and parameters in table 1 [12]. SNAP-21 is constructed so that 
it is a fully shielded device with a removable generator, making the 
handling and maintenance easier. The RTG ff s predicted reliability and 
long life fulfills the requirements for deep sea use. The maximum 
temperature of the fuel in the SNAP-21 design is 1700°F (as listed in 
table 1) while the hot end temperatures of the thermoelements is 1000°F. 
This 700°F temperature drop from the fuel pellet to the themoe lectric 
generator is due to heat conduction through the biological shield. 

A heat pipe can be introduced to the SNAP-21 design to enhance the 

^Numbers in brackets refer to bibliography 
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POWER OUTPUT 


10 watts 


EXTERIOR SURFACE MATERIAL 


Copper on Titanium 


DESIGN LIFE 


5 years 


SYSTEM OUTPUT VOLTAGE 

with power conditioner 


24 VDC 


without power conditioner 


4.7 VDC 


LOAD CURRENT 

with power conditioner 


.42 amp 


without power conditioner 


2 . 1 amp 


NET SYSTEM EFFICIENCY 
end of life 


67. 


beginning of life 


57. 


GENERATOR EFFICIENCY 


8-97. 


POWER CONDITIONER EFFICIENCY 


907. 


FUEL TYPE 


SrTi0 3 


FUEL LOADING 


33,000 curies 


FUEL TEMPERATURE MAX 


1700 °F 


T/E HOT JUNCTION TEMP, MAX 


1000 °F 


T/E COLD JUNCTION TEMP, MAX 


120 °F 


OPERATION DEPTH, MAX 


20,000 feet 


WATER TEMPERATURE 


28-41 °F 


TOTAL WEIGHT 


640 lbs 



TABLE 1 

(SNAP-21 System Parameters) 



13 



heat transfer from the heat generated in the fuel pellet to the hot plate 
of the generator. The heat pipe consists of essentially a closed evacu- 
ated pipe whose inside walls are lined with a capillary structure (called 
a wick) saturated with a volatile fluid (see figure 2). The operation of 
a heat pipe combines two principles of physics: vapor heat transfer and 

capillary action. Vapor heat transfer is used for transporting the heat 
energy from the evaporator section to the condensor section. Capillary 
action then returns the condensate back to the evaporator section. The 
working fluid absorbs the heat energy received at the evaporator section, 
transports it through the pipe, and releases it at the condensor end. 
Since evaporation and condensation take place at essentially the same 
temperature, the heat pipe operates almost isothermally [4]. For ex- 
ample, a thermal power of 11000 watts was carried 27 inches by a one 
inch diameter heat pipe with a temperature loss so small it was diffi- 
cult to measure accurately. By way of comparison, a copper block nine 
feet in diameter and weighing about 40 tons would be required to pro- 
duce the same result. [8] 
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EVAPORATOR SECTION CONDENSER SECTION 



Figure 2 

(Basic Components of the Heat Pipe) 



The objective of this thesis was to parametrically analyze the 
feasibility of adapting the heat pipe to a deep sea RTG and to compare 
the design results using the heat pipe to an existing RTG design (SNAP- 
21 ). 
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II. BASIC DESIGN OF PROPOSED RADIQISI OTOP XC THERMOELECTRIC GENERATOR 



A. COMPONENTS OF DESIGN 

Figure 3 contains a cross section of the proposed RTG system. The 
major change made to the SNAP-21 design was the implementation of a heat 
pipe to enhance the heat transfer from the heat generated in the annular 
fuel pellet to the hot plate of the annular generator. A tubular geo- 
metry was used for the heat pipe primarily because experimental data 
exists for this configuration. The annular geometry of the fuel pellet 
and thermoelectric generator was necessary to provide enough surface area 
for efficient heat transfer to and from the heat pipe. 

In order to get a meaningful comparison between designs, this 
analysis used the same materials as used in the SNAP-21 design. The bio- 
logical shield is made of depleted Uranium with 87. Molybdenum added to 
shift the crystalline phase change of Uranium from 1300°F to 1800°F. De- 
pleted Uranium is also a much better shield per weight than any other 
shielding material exhibiting similar strength. Safety requires that the 
fuel does not rupture its container under any accident. Thus, Hastelloy 
C is used as the fuel capsule due to its weldability and high tensile 
and yield strengths. Strontium-90 is the fuel used because of its rela- 
tively cheap cost and availability. The pressure vessel is made of 721 
Titanium because of its high strength per weight. The heat pipe is made 
with stainless steel containing sodium as a working fluid because of the 
expected operating temperature of 1300°F. To make the heat pipe as ef- 
fective as possible [4], a small capillary pore radius of .02 inch was 
chosen. From figure 3 of reference [4], the pressure head developed in 
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Figure 3 

(Basic features of the proposed design) 
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the sodium heat pipe was estimated to be *09 psi* These material speci- 
fications are given in table 2 [12]. 

B. IMPOSED DESIGN RESTRICTIONS 

Due to its environment* a deep sea RTG has specified conditions of 
operation* design* and fabrication not ordinarily required of such a 
device* These conditions as listed in SNAF-21 B Program quarterly re- 
port #1 [12] are : 

Temperature of surrounding water* 34.5°F + 6*5°F 

External pressure *.*. 10 ,000 psi or 22*500 feet (ocean depth) 



Humidity 100% max in salt atmosphere 

Shock an excursion producing 6 g's 



Vibration ...*an excursion producing 3 g's 
In addition to the above environmental restrictions, the desired mechani- 
cal properties of a deep sea RTG are: 

Fully shielded device 

Removable thermoelectric generator 

Life of five years 

The Atomic Energy Commission (AEG) has specified in its safety stand- 
ards the condition required for a fully shielded device. This condi- 
tion states that the external dose rate due to nuclear radiation must 
be less than 200 milli Roentgen per hour (MR/hr) [13]. 

There were no attempts made to redesign the generator thermocouple 
elements so the specifications as presented by the SNAP-21 design were 
those used in this analysis* The hot end temperature of the thermo- 
elements (T ) was held constant because the SNAP-21 design specifications 

ri 

for the thermoelements were not available for temperatures higher than 
1100°F. These thermocouple element materials are listed in table 2. 

Those factors held constant in this analysis were: the materials 

and their properties, the generator specifications (output at beginning 
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FUEL 


composition 
chemical form 
melting point 
max. operating temp 
density 

specific activity 
specific power 
power density 
specific heat 
thermal conductivity 
radiation 
half- life 


SR-90 

SrTi0 3 

1900°C 

1425°F 

3.7 gm/cc @ 70°F 

33 curies/gm 

0.21 watt/gm 

0.8 watt/cc 

0.143 BTU/gm°F 

0.0132 cal/sec-cm°C 

gamma due to bremsstrahlung 

28 years 


FUEL CAPSULE 


composition 
melting point 
thermal conductivity 
density 

tensile strength 
yield strength 


Hastelloy Alloy C 
2400°F 

10.0 BTU/hr-f t-°F @1150°F 
0.323 pounds/inch 3 Q 72°F 

121.000 psi @ 72°F 

81.000 psi @ 72°F 


HOUSING 


composition 

density 


721 Titanium ^ 

0.16 pounds/inch 


GENERATOR COLD 
FRAME 


composition 

density 

thermal conductivity 


Tellurium Copper 
0.32 pounds/inch^ 
209. BTU/hr-f t-°F 


SHIELD 


composition 

density 

thermal conductivity 


depleted Uranium 8% Moly 
.629 pounds/inch^ 

19 BTU/hr-f t-°F 


INSULATION 


composition 
thermal conductivity 


Linde Super Insulation 
0.0008 BTU/hr-f t°F 


THERMOELEMENTS. . . 


N-type 

P-type 


PbTe segmented 
(BiSb) 2 Te 3 and PbTe-SnTe 


HEAT PIPE 


wall 

fluid 

operating temperature 
capillary pore radius 
pressure head developed 


Stainless Steel 
Sodium 
1300°F 
0.02 inch 
0.09 psi 



TABLE 2 

(Materials Used in SNAP-21 Design and Proposed System) 
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of life of 13*8 watts) 9 heat pipe diameter (1 inch), fuel capsule thick- 
ness (%-inch) , pressure vessel thickness (1 inch for the cylindrical 
walls and %-inch for the hemispheres) , heat loss through shield structural 
supports and generator insulation (47 watts), and the hot end operating 
temperature of the generator (1100°F). 

Since the SNAF-21 design is an operating system and the materials 
used in the proposed design are identical with SNAF-21 specifications, 
this parametric analysis should give results comparable to the building 
of an operational system. 

C. DESIGN VARIABLES 

The parameters varied were the outer radii of the fuel, biological 
shield, and insulation (also forcing the length of the fuel pellet to 
vary) . 

Three mechanisms of the proposed RTG design (figure 3) are func- 
tions of the variable parameters. The changing of any one of the radii 
affects the safe level of the nuclear radiation at the outside of the 
pressure vessel (discussed in the Nuclear Radiation Analysis). Both the 
generator cold frame and insulation package heat transfer equations are 
basically functions of radial variations. Thus the thermoelement cold 
end operation temperature and the total system heat loss are also func- 
tions of the variable radii (these effects are clarified in section III). 
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III. ANALYSIS OF PROPOSED RTG DESIGN 



The heat transfer was analyzed using steady state, constant property, 
radial heat conduction* The nuclear radiation was treated according to 
equations given in Etherington 1 s Nuclear Engineering Handbook , SNAP-21 
specifications, and the principle of superposition of dose rates. 

A digital computer program was written to analyze this system. 

It takes into account these three effects as the outside radii of the 
biological shield, fuel pellet, and insulation are varied. Within the 
acceptable limits of the dose at the outside of the pressure vessel due 
to nuclear radiation (200 MR/hr + 25 MR/hr) , the weight of the system 
(excluding insulation and power conditioning unit) and the overall ef- 
ficiency (at beginning of life) were calculated and plotted. 

A. NUCLEAR RADIATION ANALYSIS 

The analytical solution for the dose rate at a given distance from 
a finite length cylinder of radioactive material is given in Etherington f s 
Nuclear Engineering Handbook , Chapter 7 [6]. This solution (figure 4) 
assumes a uniformly distributed source emitting radiation at only one 
energy and shielded by only one type of shielding material in the form 
of a slab. The RTG system as proposed in this design has four compli- 
cating variations. 

1) The gamma radiation from the fuel has a continuous energy 

spectrum with an E of 2.18 ev. 
max 

2) There is a series of shielding materials present -- fuel 
capsule, biological shield, pressure vessel. 

3) The fuel is an annular cylinder. 

4) The shield is an annulus circumscribing the fuel pellet. 
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FUEL 



D 





= 6 . 95 x 10 



- 5 A*-o- 



eb4> = kb4> 



where <f) - Ft e , ^tju, t c ) 

P = source strength for volume source 
yU = attenuation coefficient of shield 
/Me = attenuation coefficient of source 
t c = self absorption distance 
t = shield thickness 
- 55 —= mass absorption coefficient of air 

*10l 

E = energy of nuclear radiation 
B = buildup factor 

Figure 4 

(Model employed to calculate dose rate due to nuclear 

rad 1 at i on ) 
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The first difference was handled by breaking the energy range up 
into ten equal intervals and considering each interval as one distinct 
group at the mean energy of that group. The dose rate was then calcu- 
lated for each interval and summed over the ten energy intervals to get 
the total dose rate. The values of jj^ (attenuation coefficients) for 
Sr-90 fuel and the various shielding materials as a function of energy 
are presented in appendix A. Values of F (©,yu^u c *£c.) ^ c are plotted 
in reference [11]. 

The second difficulty was removed by replacing yut by u.ii- , 
where i refers to the different shields, and then proceeding as usual. 
Experiments made by the 3-M Company in their first quarterly report 
SNAP 21-B Program [13] indicated that this analytical model was approxi- 
mately 30% conservative. 

The third problem was overcome by assuming that the dose rate due to 
different sources can be superimposed. Using the following sketch, the dose 
rate due to A equals the dose rate due to B minus the dose rate due to C. 




Dose Rate 
Due to A 




Dose Rate 
Due to B 




C 



Dose Rate 
Due to C 
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The fourth variation was analyzed through geometrical considera- 
tions. Using the following sketch* the annular shield presents more 
shielding material to the nuclear radiation than the slab shield pro- 
posed by Etherington * s model. Thus the solution using a slab should be 




Since A ^ B, more shielding was present than that 
assumed and the proposed model was conservative. 



Etherington 1 s equation was programed on a digital computer (IBM 
360) to handle the specific case of this design. Interpolations for F 
( yu-t +yU Qm 't c ) and t£ c were done assuming linear relation- 

ships. The computer program's estimated accuracy is + 25 MR/hr (arrived 
at by checking computerized answers with hand calculated values). A 
numerical solution to Etherington ' s equation (as modified) is given in 
appendix A, as well as the computerized solution. 

B. HEAT LOSS THROUGH INSULATION 

The proposed RTG design has a series of finite length concentric 
annular cylinders restricting the heat transfer from the fuel pellet to 
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the surrounding water. In order to analyze the heat loss to the water 

it was assumed that the concentric annular cylinders were of infinite 

length. A cross section of the proposed RTG design is indicated in 

figure 5. The annular fuel pellet requires the solution of an internal 

heat generation* radial heat conduction equation. The remaining annular 

cylinders were cases of pure radial heat conduction. The analytical 

solutions for these two cases are derived in appendix B. The equations 

for the particular case of the proposed RTG were combined to find the 

heat transfer rate out, per unit length, when the inside and 

outside temperatures (T- and T.) are specified. This particular case 

1 o 

is done analytically in Appendix B. 

The adaption of this solution concerning an infinite length cylinder 
(or finite length cylinder with insulated ends) to that of a finite length 
cylinder necessitates a basic assumption. In order to make the assump- 
tion reasonable, an analysis of the ratio between the end surface area 
and the lateral surface area was made. Letting R = cylinder radius and 
L * cylinder length, the ratio of the areas becomes: 



ends 



2ttR 2 

2ttR^ 



R 



L 



lateral 

Basing the total heat transfer of the finite length cylinder on an area 
basis, the actual heat transfer (since calculated according 
to the equation given in Appendix B) is: 



l L ' 1 + stSr™ ! 

= (fL l L * r] 
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T's represent temperatures at interfaces 

kf = thermal conductivity of the fuel 

k c = thermal conductivity of the fuel capsule 

k s = thermal conductivity of the shielding material 

kj = thermal conductivity of the insulation 

Figure 5 

(Cross section of proposed design) 
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^actual ^ etween 5 and 11 watts, thus making the initial assump- 

tion of insulated ends (infinite length) close to correct. For values 
of R such that R L, the lateral surface area is much larger than 

the end surface areas. For these cases the heat transfer rate out the 
end areas would be negligible. 



C. GENERATOR COLD FRAME HEAT TRANSFER 

Tellurium Copper was the material used for the cold frame due to 
its high thermal conductivity of 209 BTU/hr-f t-°F. The basic design 
considers the heat transfer from the cold end of the thermoelements to 
the outside water as radial heat conduction with constant properties. 
The heat conduction equation for an annulus is given below [2]: 



^transferred out - W l_ 



(TI-TO) 

MRo/RO 



where RO = outer radius 
R1 s inner radius 
T1 = temperature at radius R1 
TO * temperature at radius RO 
L - length of annulus 
k = thermal conductivity 

For the assumed model shown in figure 7; 

T1 « cold end temperature of thermoelements 
R1 - 1.3 inches - radius of thermoelements 
cold end 

TO = 505 °R outside temperature of annulus 
L - 1 inch 

The temperature variation in the radial direction was approximately from 
70°F at a radius of 1.3 inches, to 45°F at the outer radius RO; over such 
a small range of temperatures the constant properties assumption is a valid 
one. 



The shape of the generator cold frame can be seen in figure 6. Since 
both the upper and lower surfaces are insulated, the heat transfer is 
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Figure 6 

(Cross section of generator cold frame) 
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Figure 7 

(Geometry assumed for heat transfer through 
generator cold frame) 
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Figure 8 

(Efficiencies vs. hot end temperature of SNAP-21 thermoelements) 




assumed to take place in the configuration shown in figure 7. The 
thermal gradients in the axial direction in the figure 7 model are 
small and the proposed model assumes them negligible. 

The efficiency ( I'j ) of a thermoelectric generator can be divided * 
into the Carnot efficiency ( IT| ) and the material efficiency (V|^) as 
indicated in figure 8 [1]. Figure 8 shows that increases with in- 

creasing temperature difference ( AT ) between the hot and cold ends of 
the thermoelements. The model used is one where V| ^ was constant and 
was thus conservative with respect to increasing AT • With 
constant, the efficiency ( Yj ) is then a function of Vj^ only, and 
with the same hot end temperature ( T^ ), is a function of 

only (or more exactly, Tl) . 

The heat conduction equation shows that the temperature Tl and the 

heat to be transferred out are linearly related. Since the efficiency 

( V) ) is a function of Tl and Q _ , , the solution is itera- 

f transferred out 

tive (i.e., pick Q . , ~ this specifies Tl, which then speci- 

9 transferred out 9 • 

fies a new Y] , which then specifies a new C) _ , t ). 

^ 9 r transferred out 7 

D. DERIVATION OF COMPUTER PROGRAM 

The subroutine HEAT was derived from the analytical model assumed 

/ 

for the heat loss through the insulation, and determines the total system 
heat loss (since all other losses were assumed constant). The subroutine 
QGEN determines the amount of fuel needed to supply 13.8 watts of power 
at beginning of life, and was derived from the analytical model assumed 
for the generator cold frame. The subroutine DOSE determines the dose 
rate in MR/hr at the surface of the pressure vessel for a given set of 
variables (developed from the Nuclear Radiation Analysis). 
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The flow chart in figure 9 shows that a given set of radii were 
picked and a total system heat loss (FUELS) was assumed. The subroutine 
QGEN then calculated the required amount of fuel needed to fulfill the 
design output (FUELG) . The total required amount of fuel was then given 
by: 

FUELT = FUELS + FUELG 

Since the actual total system heat loss was obtained from an iterative 
process, a new value for the total system heat loss (QOT) was calculated. 
QOT was then compared to the assumed value (FUELS). With a positive 
answer to the IF statement the dose rate was then calculated by the sub- 
routine DOSE, and the printable output limited to the range: 

175 MR/hr dose rate ^ 225 MR/hr 
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IV. RESULTS AND DISCUSSION 



Limiting conditions were imposed on the proposed RTG design due to 
AEC regulations specifying a safe maximum external dose rate of 200 MR/hr. 
For a given set of dimensions, the dose rate at the outer radius of the 
pressure vessel was calculated (using the computer program) with an 
assumed accuracy of + 25 MR/hr. Any information regarding total weight 
or overall efficiency of the RTG was likewise limited to a range of dose 
rate (175 MR/hr to 225 MR/hr). Most graphs have these corresponding 
ranges in dose rate plotted. Note that in the following discussion, 
weight refers to weight of major components, which is equal to total 
weight minus weight of insulation and power conditioning units. 

Figures 10 through 15 show values of weight versus fuel radius for 
constant values of insulation thickness. The information contained in 
these graphs indicates that there is a fuel radius where a minimum 
weight is reached for each value of insulation thickness. Knowing that 
the weight of various components are functions of the fuel radius and 
length (i.e., the fuel surface area), a minimum is expected. In the 
case of a solid cylinder, with a fixed volume, there is a set of dimen- 
sions where the surface area is a minimum (i.e., diameter = length). 

When values of weight versus insulation thickness at constant values 
of fuel radius are plotted (figure 16), the radius of the fuel where 
minimum weight is reached is seen to always lie between 1.3 and 1.5 
inches. When values of the fuel radius at which minimum weight is reached 
versus insulation thickness is plotted (figure 17), the general curve 
suggests using a slightly smaller fuel radius with increasing insulation 
thickness for minimum weight in the design. 
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Figure 18 plots overall efficiency at beginning of life versus fuel 
radius at constant insulation thickness and shows two aspects of the 
design: 

L Overall efficiency increases non-linearly with increasing 
insulation thickness. 

2. Maximum efficiency is reached somewhere between a fuel 
radius of 1.55 and 1.70 inches. The maximum efficiency 
is never located where the minimum weight occurs. 

Thus it was never possible to combine minimum weight with the maximum ef- 
ficiency for a specific value of insulation thickness. 

A plot of efficiency versus insulation thickness (figure 19) shows 
a non-linear increase in efficiency with insulation thickness. There- 
fore the greatest gain in efficiency resulted from the initial addition 
of .02 inches of insulation. 

In order to get some comprehensive results it was necessary to plot 
minimum weight and the efficiency at minimum weight versus insulation 
thickness. Figure 20 shows that neither the minimum weight nor the ef- 
ficiency are linear functions of insulation thickness until an insula- 
tion thickness of 1.6 inches is reached. This is recognized when the 
percentage increase in both minimum weight and overall efficiency is 
plotted versus insulation thickness as in figure 21. Thus using 1.6 
inches of insulation results in the greatest gain in efficiency at the 
lowest gain in weight. 

The SNAP-21 design has a maximum fuel temperature of 1700°F. The 
proposed RTG, through the use of the heat pipe, has lowered this tempera- 
ture to a range of from 1385°F to 1410°F. 

The SNAP-21 design operates at 5 % overall efficiency at the beginning 
of life, and has a total weight of 640 pounds. Allowing 60 pounds (as 
done in SNAP-21) for the weight of the power conditioning units and the 



34 



insulation, the design with the maximum total weight contains 613 pounds 
and has an overall efficiency at beginning of life of 5.415%. This is 
also the design with the maximum overall efficiency. The minimum total 
weight obtained for any of the proposed RTG*s was 555 pounds and had an 
overall efficiency at beginning of life of 5.37%. If the RTG was de- 
signed with 1.6 inches of insulation, the overall efficiency at begin- 
ning of life would be 5.4% with a minimum total weight of 570 pounds. 
Thus the introduction of the heat pipe has two effects, increasing the 
overall efficiency at beginning of life and decreasing the total weight. 
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Figure 10 

(Weight of major components vs. fuel radius at 1.0" insulation thickness) 
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(Weight of major components vs. fuel radius at 1.2” Insulation thickness) 
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Figure 12 

(Weight of major components vs. fuel radius at 1.4" insulation thickness) 
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Figure 13 

(Weight of major components vs. fuel radius at 1.6" insulation thickness) 
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FUEL RA0IUS (inch) 

Figure 14 

(Weight of major components vs. fuel radius at 1.8" insulation thickness) 
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Figure 15 

(Weight of major components vs. fuel radius at 2.0" insulation thickness) 
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(Overall efficiency vs. fuel radius at constant values of insulation thickness) 
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V. CONCLUSIONS 



1) The use of a heat pipe in a RTG reduces the total weight 
and increases the overall efficiency (as compared to SNAP-21), while 
maintaining reliability and long life* 

2) The design with minimum total weight has a different fuel 
radius than the design with maximum overall efficiency (at a specified 
value of insulation thickness). 

3) The use of the heat pipe decreases the maximum temperature of 
the fuel from 1700°F to approximately 14Q0 o F. 

4) In most cases, the biological shield accounts for 40% of the 
total system weight. 
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VI. RECOMMENDATIONS 



1) Perform experiments to verify the accuracy of the analytical 
solution for the dose rate due to the nuclear radiation. 

2) Parametrically analyze the system for variable operating 
temperatures of the thermoelements. 

3) Derive an analytical model for the heat transfer through the 
insulation surrounding the biological shield that doesn’t assume in- 
finite length geometry. 

4) Analyze the relationship between the higher cost but lower 
amount of shielding required for an o( emitting fuel. 

5) Consider different heat pipe geometries for a more efficient 
heat transfer system. 
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APPENDIX A 



NUCLEAR RADIATION ANALYSIS 

Considering the model whose dimensions are shown in figure 22, 
tables 3 and 4 contain all the information necessary to calculate the dose 
rate due to the annular geometry fuel pellet. The final equation for the 
dose rate is (from figure 4) [13]: 

D = 

” here ^ ' 2(Z°f to l F ^ + 

K - constant given in figure 4 
B s buildup factor 



The computer program for this specific case calculates a dose rate of 
118 MR/hr, The hand calculated dose rate is 134 MR/hr. 
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(Sample Calculations of Nuclear Radiation Dose Rate) 



From the superposition of dose rates 9 the dose rate due to A equals 
the dose rate due to B minus that due to C. Thus the dose rate in Rad/hr 
is .1714-. 06119 = .102. This is then converted into Rem/hr for the dose 
rate to be given in MR/hr (134 MR/hr) . Note that 1.31 Rem equals 1 Rad. 



APPENDIX B 



HEAT LOSS DERIVATION 



The solution to the heat transfer equation for the annular cylinder 
both with and without internal heat generation is given on the following 
page [2]. Using these solutions when considering the geometry shown in 
figure 5, there can be derived six equations with six unknowns. The 
solution of these six equations results in the calculation of 

l £ / ou t 

The six equations are" 




2TTk c 



(T, -T^ 

M'W) 



Qtt 2 ii kf 



\ 5 /OUL-i- 

I2\ - ?Tr t. (T. -TO 

ULr ^(u/r z ) 

f—\ - ? TT W CH± ~~ ) 

( l Lt " 2 s U, (T*/**) 




(Tg-Tt) 

iu ( 



( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 



The known values of k and r are then substituted to simplify the six 
equations. 



(f ), = '”< T «-TO 



( 1 ' ) 
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Annulus with no Internal Heat Generation 




k - thermal conductivity 

T " 

Q-. 2-rxk -gc T,) 

* (k( r i /{ z) 



Annulus with Internal Heat Generation 




where 



c; 




= "T-r-T? + 



Q 

4k 




+ 



4- k 



*i z ~K Z ) 

C [Hn^3 
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( 2 ') 




(fLr ' 20 -' c ' 

(fL= W.8(T r T t )i,(Vr,) 

(f L -- I IS (TV-TV) U ( •>/ r 5 ) 
(jL t = .005-02 (T 5 -T t ) (vv/r 4 ) 



O’) 

( 4 ’) 

O’) 

( 6 ’) 



The resulting six equations are then solved simultaneously. First set 
the two equations for ( S- ) , equal to one another. 

\ ft ) 

(-j) iV = \5zCX-Ti) = .012.25 Q ~ 2.0,1 c\ 



Then put in the value of C^, and solve for T^ in terms of T^ and T^. 




C^T 3 + I55T; + C, 
15-5 -f C z 



( 7 ’) 



bo^ev^e - 20.1 JUi ( . 7 5 ~ / v- 2 ') 



C? - Q 



C z 

\ 2.8 






56 



/Q \ 

Next solve for T, in terms of 7 - and T, by using equations 

3 1 j| J out 6 

4\ 5', and 6 '. 



T, = T 4 1 - C, (f )_ t 



( 8 ’) 



where - 



Ui CVO + M( r ^/r 4 ) JU (^/r 5 ) 



(aZ. & 






.OoSoZ 



'Q' 



From equations 3’ and 7* solve for in terms of ^~jp) out an< ^ 

Q 

5-5- C* 



T _j _ ( I5~5d- & )out 5 - - C-4^) - C x Cj 

^ I 



(9’) 



where 



c * = 

C 5 - 



Qtt 



( v ?) 



Q C 2 
12.8 



Then can be eliminated between equations 8 ® and 9® to find 
in terms of r^* r ^ 9 r^ s r ^ 9 and r^. 




'Q' 



out 



/ O \ _ ( 1 55 + C. y) 4 - 1 1 > b ( 1 \&) + Cg 

vT/out “ -v 15 - 5 - + c 2 



57 



APPENDIX C 



CRITICAL DIMENSIONS AND VOLUME DERIVATIONS 

Knowing the fixed outside diameter of the heat pipe (1 inch) and 
the required hot plate area of the generator (A~as specified in the 
SNAP-21 design) the length of the heat pipe needed next to the genera- 
tor to fulfill the required hot plate area is easily calculated. There 
is an additional specification of 1/16 of an inch gap between the heat 
pipe and the hot plate for heat transfer by conduction and radiation 
(at a temperature drop of 200°F). Figure 23 shows the geometry and the 
required hot plate length. 

From the figures 24, 25, 26 and 27 the respective volumes of the 
generator cold frame, pressure vessel, shield, and fuel capsule can be 
easily verified. Note that in each case these volumes are functions of 
the variable radii (RF-radius of the fuel, RS-radius of the shield, 
Rl-radius of the insulation). 
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total hot plate area = A = 6.6 in^ 

A s JX. ° 2 L D = 1 1 In 

L » 6.7 in 

Figure 23 

(Length dimension of annular generator) 
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Figure 24 

(Volume derivation of annular generator) 



VOLUME 




(4/3) ((Rl + l. ) 3 -(R| + . 5)3)^] 



Figure 25 

(Volume derivation of pressure vessel) 
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where 1 tRS^( LF+ . 5+2 ( RS-A ) ) * volume of solid 
cyl Inder 

TT( RS-A ) ( 1/ 2 ) 2 = volume cut out for 
passage of heat pipe 

TTA (LF+.5) = volume cut out for fuel and 
fuel capsule 

TT. 86A(RS-A)2 = volume cut out for the 
rounded corners (i.e. see below) 




( 1-Tty 4 ) (RS-A ) 2 

volume * 2TT (2A) (area) 

* (4-7T )7TA(RS-A) 2 

= . 86A ( RS-A ) 2 



Figure 26 

(Volume derivation of shield) 
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Figure 27 

(Volume derivation of fuel capsule) 
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